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We use scanning tunneling microscopy to determine the surface structure and dopant distribution 
in Pr cc Cai_ x Fe2As2, the highest-T c member of the 122 family of iron-based superconductors. We 
identify the cleaved surface termination by mapping the local tunneling barrier height, related to the 
work function. We image the individual Pr dopants responsible for superconductivity, and show that 
they do not cluster, but in fact repel each other at short length scales. We therefore suggest that the 
low volume fraction high-T c superconducting phase is unlikely to originate from Pr inhomogeneity. 

PACS numbers: 68.37.Ef, 74.55. +v, 74.70.Xa, 74.62.Dh 



The recent discovery of high-T c superconductivity in 
Fe-based materials [1] has rejuvenated worldwide efforts 
to understand and predict new superconductors. Fe- 
based superconductors (Fe-SCs) share several character- 
istics with their cuprate cousins. First, their chemical 
structure consists of Fe-based superconducting planes 
separated by buffer layers. Second, superconductivity 
usually arises by electron or hole doping a stoichiometric 
antiferromagnetic parent compound [2 . Third, Fe-SCs 
show nanoscale electronic inhomogeneity [3j 0], necessi- 
tating the use of local tools for complete understanding. 

The first generation of the AFe2As2 (122) family of 
Fe-SCs seemed similar to cuprates, because hole dop- 
ing this family resulted in a much higher maximum 
T c (38 K in K a; Bai_ :c Fe2As2 [5]) than electron doping 
(25 K in Ba(Fei_ :E Co :c )2As2 [6 ). It was pointed out 
that the highest T c of all Fe-SCs arose in an electron- 
doped member of a different structural family (57 K in 
Smi-^La^Oi-^F^FeAs [7 ). This suggested that the T c 
in electron-doped 74(Fei_ :E Co :E )2As2 might be artificially 
suppressed by the Co dopants, which play a dual role of 
providing the charge carriers necessary for superconduc- 
tivity, but also introducing damaging disorder directly 
into the crucial Fe layer. It was hypothesized that elec- 
tron doping into the buffer layer between the crucial Fe 
layers could provide the charge carriers without the disor- 
der, and perhaps allow a higher T c in the electron-doped 
122 materials. Indeed, this was found to be the case in 
the rare-earth-doped Cal22 family [8j [9], with T c reach- 
ing 49 K in Pr x Cai_ :c Fe2As2. However, the high T c ap- 
peared in only ~ 10% of the volume, while the bulk of 
the material showed T c ~ 10 — 20 K. 

Saha et al. performed a thorough search for the ori- 
gin of the low volume fraction high-T c phase, using bulk 
experimental probes. First, the high T c was found to be 
impervious to etching or oxidation, arguing against sur- 
face superconductivity. Second, high-T c resistive tran- 
sitions were never observed for dopant concentrations 
below those necessary to suppress the parent antiferro- 
magnetic phase, arguing against any randomly occurring 



impurities or contaminant phases as the origin. Fur- 
thermore, no such contaminant phases were observed in 
over 20 samples examined by x-ray diffraction. Third, 
the high T c was unaffected by the global structural col- 
lapse phase transition (the abrupt ~ 10% shrinkage of 
the c-axis lattice constant that occurs in the Cal22 fam- 
ily under external or chemical pressure), arguing against 
any relationship to the collapsed phase or to interfaces 
between collapsed and non-collapsed phases. Finally, it 
was noted that aliovalently-doped CaFe2(Asi_ ;E P ;E )2 also 
shows the structural collapse but no high-T c volume frac- 
tion [10] . Saha et al. therefore concluded that the charge 
doping is an essential ingredient to the high-T c phase, 
and speculated on the low volume fraction high-T c phase 
"as having a localized nature tied to the low percentage 
of rare earth substitution." 

Given the challenges in identifying the origin of the low 
volume fraction high-T c phase from bulk experiments, a 
local probe is naturally required. Here we use scanning 
tunneling microscopy (STM) to investigate two possible 
sources of electronic inhomogeneity in Pr iC Cai_ iC Fe2As2: 
the surface and the dopants. We provide the first defini- 
tive identification of the cleaved surface termination and 
the first image of all individual dopants in the Cal22 sys- 
tem. Our results suggest that dopant inhomogeneity is 
unlikely to be responsible for the low volume fraction of 
high-T c superconductivity in Pr x Cai_ ;E Fe2As2. 

Single crystals of Pr x Cai_ :E Fe2As2 are grown via self- 
flux with measured x = 10.5% and resistive T c = 43.2 K. 
The crystals are handled exclusively in Ar environment, 
cleaved in ultra- high vacuum at cryogenic temperature, 
and immediately inserted into the STM head where 
they are imaged with a mechanically sharpened Ptlr tip, 
cleaned by field emission on Au. The first challenge in 
STM imaging of any new material is to identify the sur- 
face structure and evaluate to what extent it is represen- 
tative of the bulk. The surface structure of the AFe2As2 
system has been particularly controversial [TT] , Due to 
the stronger bonding within the FeAs layer [Fig. 1(a)], 
the FeAs layer is expected to remain intact, leaving half 
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FIG. 1. (color online) Surface morphologies of cold- 
cleaved Pro.io5Cao.895Fe2As2. (a) Crystal structure of 
Pr x Cai_ x Fe2As2. Topographs of (b) 2 x 1 surface structure 
(250 pA, +300 mV, 7 K) (c) disordered, "web-like" surface 
structure (15 pA, 100 mV, 7 K), (d) smooth transition be- 
tween 2x1 and "web-like" structures (20pA, -lOOmV, 7K), 
and (e) 1 x 1 square lattice with ~ 4 A lattice constant (5 pA, 
50 mV, 25 K). Inset in (b) shows an average 30 x 15 pixel, 2x1 
supercell [17] tiled 2x4 times. Inset in (e) shows an enlarged, 
4 nm topograph of 1 x 1 square lattice acquired at 50 mV and 
50 pA. (Due to an external noise source present during the 
acquisition of the data in panel (e), the images in (e) have 
been filtered to remove all spurious spatial frequencies higher 
than the 4 A periodicity.) 



TABLE I. Work functions for several elements [18] . 
Atom Fe As Ca Pr Sr Ba Ail Pt Ir 
<p (eV) 4.65 3.75 2.71 2.7 2.76 2.35 5.32 5.40 5.6 

a complete A layer on each surface [12]. This preserva- 
tion of charge neutrality is a necessary (but not sufficient) 
condition for the surface to be representative of the bulk. 
However, a number of experiments have claimed that the 
cleaved surface is As-terminated in the Bal22 [13] . Srl22 
[T4] . and Cal22 [15] systems. 

We encounter three different surface morphologies in 
our STM topographs of Pro. 105 Cao. 895^2 As2- The ma- 
jority of the observed sample surface displays a 2 x 1 
structure [Fig. 1(b)] frequently observed in other STM 
studies of 122 materials [16]. We occasionally observe 
a disordered, "web-like" structure [Fig. 1(c)], which 
smoothly merges with the 2x1 structure [Fig. 1(d)]. 
The third type of surface, observed rarely, shows a 1 x 1 
square lattice with ~ 4 A periodicity [Fig. 1(e)]. 

We map the tunneling barrier height to identify these 
surfaces. The tunneling current / is expected to decay 
exponentially with the tip-sample separation z as 

/ 8m e ¥ 
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where is the local barrier height (LBH) , approximately 
equal to the average of the tip and sample work func- 
tions [19 . Due to the large variation in the work func- 
tions of the elemental components of 122 Fe-SCs (Table 
|l|, a simple measurement of the LBH should in princi- 
ple be sufficient to determine the surface termination. 
However, the LBH is sensitive not only to the elemental 
composition of the tip, but also the geometric configu- 
ration of the tip's terminal atoms, and the tip-sample 
angle (which reduces the LBH by cos 2 where is the 
deviation between the sample surface perpendicular and 
the z direction of tip piezo motion [20]). Moreover, the 
LBH depends on the sample topography through two op- 
posing mechanisms. On the one hand, protruding atoms 
or clusters may deform easily, stretching out as the tip 
is retracted, decreasing the effective rate at which the 
tip-sample distance decreases, and thus suppressing the 
measured LBH above the protrusion [21]. On the other 
hand, the topographic corrugation appears smoothed out 
at distances far from the surface; this implies that the 
wave function must decay faster above a protrusion than 
a depression, thus enhancing the measured LBH above 
a protrusion [22 . Without accounting for these factors, 
previous studies found the LBH on the 2x1 surface of 
BaFe2 As2 to be much lower than the expected work func- 
tions for either Ba or As [23 . 

In contrast, the comparison of LBH measurements 
with the same tip (i.e. the same microscopic configura- 
tion of terminating atoms) across different flat regions of 
the same cleaved surface (i.e. the same tip-sample angle) 
can yield a robust measure of relative work functions, and 
can be utilized for element identification in cases where 
the sample consists of two different surfaces [24]. Here, 
we directly compare LBH values measured with the same 
STM tip across the different morphologies of Fig. 1 on the 
same cleaved sample. 

In order to extract the LBH, at each point (x, y) in a 
field of view (FOV) a feedback loop adjusts zo(x,y) to 
maintain / = 100 pA at V^ et = —100 mV; the feedback 
loop is then discontinued and the current I(z) is mea- 
sured as the tip is retracted continuously from zo. Fig- 
ure 2 shows simultaneous topographs and work function 
maps for the lxl square lattice across a step edge [Figs. 
2(a,c)] and the 2x1 structure in a nearby fiat area [Figs. 
2(b,d)]. Figure 2(e) shows two sets of representative I(z) 
curves from the square regions in Figs. 2(a,b). Despite 
the small spread within each set of I(z) curves, the two 
sets are clearly distinct from one another. After correct- 
ing for the surface slope, we find the average LBH values 
are $ lxl = 4.50±0.42 eV and $ 2x i = 3.57±0.34 eV. In 
comparison with Table [TJ and after taking the tip work 
function into account, these values suggest that the lxl 
surface is a complete As layer, while the 2x1 surface is 
a half Ca layer. 

To quantify this identification, we note that charge re- 
distribution at a surface may give rise to an additional 



dipole barrier, which can increase or decrease the mea- 
sured LBH according to the dipole orientation [25]. We 
would therefore expect 

^As + ^d,As + ^tip -r _ ^Ca + ^d,Ca + <£tip 

^ixi — ^ ' ^2xi — ^ 

(2) 

where cp\ s and cpc a are the element work functions, 
and £?d,As and £?d,Ca are the additional energies for 
an electron to escape the dipole layers at the As- and 
Ca-terminated surfaces. Assuming (/? tip < (/?p t i r [26] . 
Eq. ([2| yields values of E^As? ^d,Ca, and their difference 
£?d,As — ^d,Ca = 0-82 eV that are all of the correct magni- 
tude for such dipole layers [25]. The sign of the difference, 
which indicates that it is harder to remove an electron 
from the dipole barrier of the As surface than that of the 
half-Ca surface, is physically justified because the half-Ca 
surface is nonpolar, whereas the As surface is deficient of 
electrons from the stripped Ca, and thus more electroneg- 
ative. The inferred electron-deficiency of this As surface 
is consistent with the failure to observe even proximity- 
induced superconductivity on the As-terminated surface 
of the related Sro.75Ko.25Fe2As2 [27] . 

We also note the reduced LBH along the step edge in 
Fig. 2(c). Such LBH reduction at step edges has been 
seen in other systems, and may be attributed to a com- 
bination of two mechanisms [24] . First, the step edge 
is effectively an angled surface, so the LBH is reduced 
by cos 2 6. Second, the Smoluchowski smoothing of the 
electron wave functions along the step edge results in an 
additional dipole moment which reduces the LBH [28 . 

To further support the identification of the 2x1 sur- 
face, we show a high-resolution map of the intra-unit cell 
structure to rule out the possibility of "hidden" surface 
atoms. We correct for small piezoelectric and thermal 
drift by placing the Ca/Pr atoms of Fig. 1(b) on a per- 
fect lattice [29] . We then use the whole FOV to create 
the high-resolution average 2x1 supercell in the inset 
to Fig. 1(b) [T7] , We do not observe atom dimerization 
(as seen in Cao.83Lao.i7Fe2As2 [IS] and Sri_ cc K ;E Fe2As2 
[HJ [30]), but rather a single row of atoms, similar to the 
CaFe2As2 parent compound [51] . 

For completeness, we investigate the nature of the 
"web-like" surface. Because it merges smoothly into the 
2x1 surface without any sharp transitions that could 
be attributed to step edges [Fig. 1(d)], the "web-like" 
surface is also likely a reconstruction of the Ca layer. Si- 
multaneous topograph and LBH map of the "web-like" 
surface are shown in Figs. 3(a,c), with analogous maps 
for the 2x1 surface, acquired with the same tip for direct 
comparison, shown in Figs. 3(b,d). Bright spots in the to- 
pograph of Fig. 3(a) exhibit much higher LBH than sur- 
rounding areas. Anomalously high LBH at these atomic- 
scale features reinforces the importance of the geometric 
effects of protrusions on the measured LBH, and the ne- 
cessity of fiat atomic planes in order to extract a reliable 
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FIG. 2. (color online) LBH comparison between lxl and 
2x1 surface structures. Topographs acquired at 25 K of (a) 
1 x 1 (4 A x 4 A) square lattice appearing on both sides of a 
step edge and (b) 2 x 1 (8 A x 4 A) lattice. Simultaneously 
acquired LBH maps are shown in (c) and (d). Approximately 
26% larger LBH is observed in a clean, flat area of the lxl 
surface (red box in (a)) than in a clean flat area of the 2x1 
surface (green box in (b)). Furthermore, sparse topographic 
protrusions on the lxl surface (e.g. marked by yellow ar- 
row) show a lower LBH close to that of the flat 2x1 surface, 
suggesting that they are scattered remaining Ca or Pr atoms. 
Both datasets were acquired at 7 S et = 100 pA and \4 e t=-100 
mV. (e) Representative sets of I(z) curves from square re- 
gions in (a) and (b) are shown as thin red and green lines 
respectively. Darker red and green lines represent linear fits 
to the average I(z) curves from boxes in (a) and (b). 



LBH comparison. Further conclusions from Fig. 3 would 
require a quantitative understanding of these geometric 
protrusion effects. We reiterate that our identification of 
the lxl surface as a complete As layer and the 2x1 
surface as a half-Ca layer is robustly drawn from the flat 
surfaces in Fig. 2. 

Since ipc& and cpp r differ by less than 1% (Table [I]), 
LBH mapping cannot be used to identify Pr atoms in 
the Ca surface layer. However, STM can image dopants 
as atomic-scale features in the differential conductance 



b 


; h>- d 


















\- • 

4nm 




i — 


,> 






— 1 High 2eV^H 


16.5 eV 



FIG. 3. (color online) LBH comparison between "web-like" 
and 2x1 surface structures. Topographs acquired at 7 K of 
(a) "web- like" surface and (b) 2x1 surface. Simultaneously 
acquired LBH maps are shown in (c) and (d). Anomalously 
high LBH values are measured over the bright, atomic-scale 
features in (a). Both datasets were acquired at 7 S et=105 pA 
and Vset—lOO mV. The z calibration used here was obtained 
by assuming that the average LBH for the 2x1 surface in the 
boxed region of (b) here is the same as that in Fig. 2(d). 



dl/dV, which is proportional to the local density of states 
[32] . Substituting Pr 3+ for Ca 2+ creates a localized pos- 
itive charge, so the impurity state is expected above the 
Fermi level. We therefore search for Pr dopants in dl/dV 
images at high bias. Figure 4(a) shows a dl/dV image 
obtained simultaneously with the topograph in Fig. 1(b) 
at +300 mV, revealing a set of bright, atomic-scale fea- 
tures. These features, which start to appear in dl/dV 
at biases higher than +70 mV, comprise ~10.4% of the 
total number of visible atoms in this FOV [33], match- 
ing the macroscopically measured x = 10.5%. The ap- 
pearance of the expected concentration of Pr 3+ dopants 
on the surface confirms the half-Ca termination. Al- 
though a subset of Co dopants were previously imaged 
in Ca(Fei_ :c Co :E )2As2 [31], this is the first time that all 
dopants have been imaged in a Cal22 system. 

Because we have imaged all dopants, we can investi- 
gate the possibility of clustering, which was suggested 
as the origin of the inhomogeneous high-T c phase [8]. 
We compute a "radial distribution ratio" (RDR) by his- 
togramming all observed Pr-Pr distances within a FOV, 
then dividing this observed histogram by an average his- 
togram of 1000 simulated random dopant distributions at 
the same concentration. The RDR in Fig. 4(b) shows no 
clustering, and in fact slight repulsion of the Pr dopants 



FIG. 4. (color online) Map of Pr dopants, (a) dl/dV image at 
+300 mV showing Pr dopants as bright atomic-scale features 
from the same FOV as Fig. 1(b). (b) Radial distribution 
ratios for two sets of Pr dopants. Full squares represent the 
distribution of Pr dopants shown in (a), while open squares 
represent a different dataset used to confirm the conclusions. 
Inset shows a 2.5 nm x 2.5 nm region in which surface Ca 
positions (white dots) and Pr dopants (yellow dots) have been 
marked, demonstrating our ability to resolve individual Pr 
dopants even at adjacent Ca sites. 



at short distances, possibly due to their like charges. The 
repulsion is not an artifact of poor dopant identifica- 
tion, as illustrated by clear detection of two adjacent Pr 
dopants in the inset to Fig. 4(b). The lack of dopant clus- 
tering in Pr ;E Cai_ ;E Fe2As2 contrasts with the Se dopants 
in FeTei-^Se^ that are prone to forming patches of ~1 
nm size [34]. This contrast may arise from the ~10% size 
mismatch of Se (198 pm) and Te (221 pm) vs. the simi- 
lar sizes of Ca (126 pm) and Pr (126.6 pm) [27j[35]. Our 
observation of the expected number of Pr dopants, more 
homogeneously distributed than would be expected for a 
random distribution, suggests that dopant clustering is 
unlikely to be responsible for the small volume fraction 
high-T c superconducting state. 

In conclusion, our STM images of the highest-T c mem- 
ber of the 122 family, Pr ;E Cai_ ;r Fe2As2, have addressed 
its surface structure and dopant distribution, with bear- 
ing on the high-T c volume fraction. First, we used LBH 
mapping to identify the 2x1 surface as a half-Ca termi- 
nation, and the lxl surface as an As termination. This 
LBH mapping method could be used to resolve debated 
cleaved surface terminations in a wide variety of mate- 
rials, such as other Fe-SCs [36] or a number of heavy 
fermion materials which have recently been cleaved and 
imaged by STM for the first time [37j [38]. Second, we 
imaged the Pr dopants responsible for superconductiv- 
ity, the first time that all dopants have been identified in 
Cal22. We demonstrated that Pr dopants do not cluster, 
and in fact show a slight tendency to repel each other at 
very short length scales. Our findings suggest that Pr 
inhomogeneity is unlikely to be the source of the high-T c 
volume fraction, in contrast to previous speculation [8]. 
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